UNCLASSIFIED 


r 


Reproduced. 
Ity.  the. 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  sped* 
flcations  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


PRESSURE  DISTRIBUTION 

ON  A  WEDGE  ACCELERATED  IMPULSIVELY 
AT  A  SUPERSONIC  MACH  NUMBER 

DECEMBER  28, 1962  DOUGLAS  REPORT  SM-42649 


MISSILE  8c  SPACE  SYSTEMS  DIVISION 
DOUGLAS  AIRCRAFT  COMPANY,  INC. 

SANTA  MONICA  CALIFORNIA 


‘w,..  \ 


f 


PRESSURE  DISTRIBUTION 

ON  A  WEDGE  ACCELERATED  IMPULSIVELY 
AT  A  SUPERSONIC  MACH  NUMBER 

DECEMBER  28, 1962 
DOUGLAS  REPORT  SM-42649 


Approved  by: 


Chief,  Micelle  Aefe/Thennodynomice  Section 


Prepared  by: 
A.I.  Ormsbee 
Consultant, 

Missile  Aero/Thermodynamics  Section 


PREPARED  UNDER  THE  SPONSORSHIP  OP  THE 
DOUGLAS  AIRCRAFT  COMPANY  INDEPENDENT 
RESEARCH  AND  DEVELOPMENT  PROGRAM 
ACCOUNT  NUMRER  ISO  10076 


MISSILE  A  SPACE  SYSTEMS  DIVISION 
OOUOLAS  AIRCRAFT  COMPANY.  INC. 


MISSILE  SYSTEMS  ENGINEERING 


fjf’l/ft 


^  SfA'  I  9  «  •  -  ». 

«>'  <  r>  M  r  *  f  .’  V. 

*.  A  N  *A  i  ,*s.  i .  A  i  A,i»  '  ** 

&-&oO~HjAT-AM-c, 

Raljscrt :  Ti aa—lttol  of  ftsofl —  Alroroft  Caniiy,  lac. ,  JUssorcfc 

Reports  MJifcU03 ,  ut  9MIM9 

To:  laodtooriors 

Mnoad  Carrie—  RwtolMJ.  lafoxaotloa  Afoaey 
Mr  Mm  Systo—  Comma* 

Ibltod  Mtw  Mr  fore* 

M-nut—  MU  Matt— 
irllja|te>  12,  Tlnlala 

Afctaartlooj  TZFR 

Josaffc  Mol 

(kl«f,  AtMUltM  kfMh 
iMMMl  hrUMIlfll  MrlSlao 

1.  It  id  rafaaadsi  that  tM  fsU«vU|  aaala— lflad  reports 

to  UHsi  so  asp  nsyisli  pMl«stl«  UaU  or  Irion  Mstrlbotod 
kgr  MO,  aw*  a*  Mo  Mriiml  Metro* t  MLUUi.  MatrlWtlsa  of 
Mism  ispota  to  oar  iataroatod  par—  or  ar^olaotl—  la  tool  rod. 


ABSTRACT 


I 


The  nonsteady  flow  equations  associated  with  an  Impulsively  accelerated 
supersonic  wedge  were  transformed  Into  £  time -independent  system.  Oils 
system  was  then  linearized  vtth  respect  to  the  velocity  change.  It  was 
shown  that  for  certain  reglonrspf  tjfe  flow  field  a  velocity  potential 
could  he  defined,  and  a  formal  sfeftes  solution  for  the  potential  In  these 
regions  was  found.  It  was  also^hotm  that  the  pressure  perturbation 
satisfies  Laplace's  equation,  Jin  a  suitably  transformed  space,  for  both 
rotational  and  lrrotational  /eglons.  Anally,  for  the  special  case  of 
a  small  wedge  angle,  the  coefficients  In  the  velocity  potential  series 
expansion  were  determine?  and  a  closed-form \gcpresaion  for  the  pressure 
field  was  obtained. 
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NOMENCLATURE 


Symbol 

o  Acoustic  speed 

Specific  heat  at  constant  pressure 
Direction  normal  to  a  streamline 

p 

Pressure 


♦  Time 

U,  Uj  Velocity  components  in  X, ,  X2  plane 

V  Initial  wedge  velocity 

AV  Increment  in  wedge  velocity 

•'i  Dimensionless  perturbation  velocity 

w 

l  Reference  Mach  nunber 

X|,X2  Physical  space  coordinates 

a  Wedge  angle 

0  Shock  angle  relative  to  wedge  velocity 

?  Ratio  of  specific  heats 


KCMBHCLATURB 


Symbol 


ton 


.1 


*2 

*1 


A 


Dimensionless  perturbation  pressure 


*1 


ii 

•l* 


f  Density 

9  Dimensionless  perturbation  density 

^  Velocity  potential 


1.  INTRODUCTION 


Current  design  philosophy  in  missiles  and  space  craft  has  generated  an 
interest  in  the  effects  of  higi  longitudinal  accelerations.  The  effects 
of  lateral  accelerations  have  been  studied  rather  extensively  (Refs.  1, 

2,  3),  particularly  within  the  framework  of  linearized  theory.  However, 
little  has  been  done  to  determine  longitudinal  effects. 

Dimensional  considerations  and  order  of  magnitude  studies  indicate  (Ref. 
l)  that  small  values  of  the  parameter  ql/V8,  where  q  is  acceleration,  I 
is  characteristic  length,  and  V  is  velocity,  should  correspond  to  small 
acceleration  effects.  Since  boost  accelerations  of  the  order  of  several 
hundred  g  are  contemplated  for  existing  designs,  and  since  these  high 
accelerations  will  prevail  through  the  low  supersonic  end  of  the  speed 
range,  a  quantitative  measure  of  the  aerodynamic  effects  is  desirable. 

Uhls  study  attempts  to  provide  insight  into  the  flow  fields  associated 
with  accelerating  bodies  by  considering  the  case  of  a  two-dimensional 
wedge  which  undergoes  an  instantaneous  but  small  velocity  change  from 
one  supersonic  speed  to  another. 

Following  Spitzer  (Ref.  4),  the  nonsteady  flow  equations  for  this  problem 
are  transformed  into  a  time -Independent  system.  A  linearization  of  the 
resulting  equations  is  effected  by  ignoring  terms  of  second  and  higher 
orders  in  the  velocity  change,  AV,  considering  the  flew  field  as  a  per¬ 
turbation  fran  the  steady  attached  shock  wedge -flow  case.  The  resulting 
linear  system  is  then  examined  from  the  standpoint  of  possible  solution, 
including  the  derivation  of  a  potential  equation  for  the  case  of  irrota- 
tional  flow.  It  is  shown  that  for  the  general  case  the  pressure  is  a 


solution  to  Iaplace's  aquation  In  a  suitably  transformed  space,  and  solu¬ 
tions  for  the  special  case  of  a  small  wedge  angle  are  obtained. 

2.  WAVE  GEOMETRY 

Consider  a  wedge  (Fig.  l)  with  initial  velocity  V  which  Is  increased 
Instantaneously  at  t-0  to  V*AV  t  with 

•A*L«1 

V 

We  choose  a  coordinate  system  ( Xj,  Xj )  which  is  translating  at  the  final 
speed  of  the  wedge  with  the  origin  coincident  with  the  wedge  vertex  for 
»*0  . 

The  wedge  velocity  is  assuned  sufficiently  large  that  a  plane  shock  wave 
is  attached  at  the  wedge  vertex.  Since  neither  the  wedge  nor  the  gas 
have  a  characteristic  length,  one  must  suppose  that  for  any  time  t  a 
length  defined,  for  example,  by  Vt  will  serve  to  define  the  scale  of  the 
flow.  That  is,  the  application  of  an  impulsive  acceleration  to  the  wedge 
generates  a  wave  pattern  which  grcera  uniformly  with  time.  The  diameter 
of  this  wave  pattern  becomes  clear  on  further  investigation.  For  a  point 
on  the  body  remote  from  the  vertex,  the  impulsive  acceleration  generates 
a  wave  propagating  in  a  direction  normal  to  the  wedge  surface  which  in¬ 
duces  a  normal  velocity  AV  sin  a  in  the  gas.  Far  AV  small,  the  propaga¬ 
tion  speed,  relative  to  the  gas,  of  this  wave  will  be  s(  ,  the  acoustic 
velocity  behind  the  initial  shock  wave .  Simultaneously,  a  cylindrical 
wave  will  originate  at  the  vertex  at  t-0  .  This  wave,  for  smell  AV, 
will  propagate  radially  from  its  center  at  speed  »| ,  and  its  center  will 
convect  along  the  wedge  surface  at  the  speed  of  the  gas.  A  new  shock 
wave  will  be  formed  at  the  vertex,  corresponding  to  steady-state  flow  at 
the  new  wedge  velocity.  The  wave  geometry  for  a  particular  time  la  shewn 
in  Fig.  2 .  The  dashed  portion  BC  of  the  shock  fails  inside  the  circle  de¬ 
fining  the  cylindrical  wave;  its  shape  is  not  known.  If  a  is  nsstmmd 
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sufficiently  small,  then  OB  and  CC'  will  he  tangent  to  the  circle.*  The 
geometry  shown  is  for  a  relatively  low  Mach  number  so  that  DD*  is  tangent 
to  the  circle.  For  sufficiently  high  Mach  numbers,  C  and  D  coalesce  and 
this  point  moves  down  the  right  hand  side  of  the  circle  toward  E  as  the 
Mach  number  increases. 

3.  TRANSFORMATION  OF  THE  EQUATIONS 


The  momentum,  continuity,  and  energy  equations  can  be  expressed  for  this 
problem  as 


Jui  ♦  u,  _fii _ i  it 

9\  p  (?X| 


j£_.  iiiSL.o 

at  a  *| 


i,j«  1.2 


(i) 


(2) 


D  /  o2 
d»  'y-i 


j  ' » 


> 


(3) 


where  it  is  assumed  that  the  fluid  is  a  perfect  gas. 
The  previous  discussion  suggests  the  transformation 


(4) 


If  the  flow  variables  are  functions  only  of  ,  Eqs.  (l),  (2),  and  (3) 
become 


* 

It  was  pointed  out  to  the 
angle  B  FCis  proportional 
restriction  on  this  angle 


author  by  Professor  Nicholas  Rott  that  the 
to  a  H  for  small  a,  placing  a  more  severe 
than  would  at  first  seem  necessary. 


<Uj 
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(5) 


li  ♦  _L  It.  - 

*(j  P  9t\ 


d  (p  U, ) 
*<\ 


(6) 


(7) 


We  now  define 


v 


I 


(8) 


X  ■ 


p-pl 

ypi 


(9) 


tm. 

pl 


(10) 


*i "  *1  -  wj 


(11) 


where  Wj  is  the  (constant)  vector  Mach  number  behind  the  initial  shock 
wave  and  ,  Pi  are  the  pressure  and  density,  respectively,  behind  the 
initial  shock. 


In  the  %  ,  9  a  coordinate  system  the  cylindrical  wave  ABCDE  maps  into  half 
of  the  unit  circle  with  F  mapping  into  the  origin.  The  waves  and  the 
wedge  are  stationary  in  this  space  (Fig.  3)* 
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We  formally  assume  that  the  restriction  of  small  AV/V  implies  small  values 
for  vj  ,  a ,  andX  .  Inserting  Eqs.  (8),  (9),  (10),  and  (ll)  into  Eqs.  (5), 
(6),  and  (7)  and  retaining  only  lowest  order  terms  provides 


*j±L 

-  J±- 

J  9n j 

9  >»  | 

’J-r- 

-In.- 

*  9  nj 

9nj 

[(y-n 

9  9  ijj  9lj 

Eqs.  (12)  and  (13)  may  he  substituted  into  (l4)  to  provide 


(12) 


(13) 


(14) 


•>!  ’J 


9  vi 


’i 


-  in.  -0 

9  i»j 


(15) 


We  consider  Eqs.  (l2),  (13),  and  (15)  as  our  system,  noting  that  Eqs.  (12) 
and  (15 )  may  be  considered  Independently  of  (l3 )  • 


’ll 


FIGURE  3 


4.  SOLUTION  OF  SHE  EQUATIONS 


A  nmber  of  manipulations  can  be  performed  on  this  system.  It  is  shown 
in  Ref.  4  that  the  characteristic  directions  for  the  systems  are  real  out 
side  the  unit  circle  and  imaginary  inside.  Since  the  solution  outside 
the  unit  circle  can  be  obtained  from  ordinary  shock  wave  considerations, 
knowledge  of  the  characteristic  system  does  not  seem  to  have  much  real 
significance  except  to  define  formally  the  unit  circle  as  the  boundary 
between  the  elliptic  and  hyperbolic  regions  of  the  plane.  Examining  Pig. 
3,  we  see  that  the  "streamlines"  throu^a  points  B  and  C  must  terminate  at 
F,  defining  a  triangular  shaped  region  outside  of  which  the  fluid  parti¬ 
cles  are  acted  upon  by  a  plane  shock  only.  The  flow  in  regions  ABF  and 
inFCE  will  thus  be  irrotational.  If  the  interaction  BC  is  a  "smooth" 
one,  then  the  entropy  gradient  d»/dn  and,  hence,  by  Crocco's  Theorem,  the 
vorticity,  will  be  of  order  A  V  in  BCF . 


If  a  velocity  potential  is  postulated 


Vi  m  j.f  ■ 


) 


it  can  be  shown  that 


U6> 


(,»  -Dili  ♦  a„  »,  £±~  ♦  (I,1  -  D-ii-  -0  (17) 

Solutions  to  this  equation,  using  separation  of  variables,  can  be  expres¬ 
sed  as 


*•  f  [\»<a,r)  ♦  B„»  (-n,r)][c;eo#nrf  *  D'„  •!«»  «*] 


(18) 
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where 


(19) 


r  -  <  *f  i2  ♦  nj  )h 

*-  ton  )  (20) 

X  ’ll 

and  f  ( n*  f)  is  expressed  in  terms  of  the  hypergeanetric  function 

Jl-.JLzi;  «+  1  it2) 

fUr)  -r"  - 2 - -  (21) 

F^i.JL^JL  ;  «♦  1;  1) 

Die  presence  of  an  infinity  in  f(-n,r)  at  r  -  0  suggests  that 

B„  -  0  (22) 

i 

We  may  also  choose,  with  no  loss  in  generality, 

1 

!  >  (23) 

1 

giving 

! 

4  -  X  f(n,  r)  (Cj,  om  fi  #+  D'„  sin  n  0)  (24) 

I 

|  These  solutions  apply  in  regions  ABF  and  FCE. 

5.  THE  PRESSURE  EQUATION 

With  same  manipulation,  the  velocities  may  be  eliminated  from  Eqs .  (12) 

|  and  (15),  providing 

I  1 


8 


j.  i»  «  d2 A 
+  2>I i  >?2  r - - — 


+  2,a.iL-  -0  (25) 
*| h 


(v?  +  ])Ji± 


♦  <,,*•►  i)iii_  ♦  2,,  JJL. 

dif22  3 1?, 


or,  in  polar  coordinates, 


r2  (r2  +  1) 


+  r(2r2  +  1)il-  + 
dr 


a2*  .  o 
a*2 


(26) 


The  transformation 

^  .vfTT2 -i 

r 


transfonns  Bq..  (26)  into  Laplace’s  equation 

t1  -iii.  ♦  r  J±-  *  £±-  -  o 
rr*  »t  «> 

with  general  solution 

X  -  +  9 


(27) 


(28) 


(29) 


6.  BOUNDARY  CONDITIONS 


The  flow  properties  in  the  region  exterior  to  the  unit  circle  are  known, 
so  that  the  values  of  pressure  and  velocity  perturbations  on  segments  AB , 
CD,  OE,  as  well  as  the  velocity  and  normal  derivative  of  the  pressure  at 
the  wedge  surface  are  known. 


On  A B 


v\ 


AV  tin  a  tin  /92 
eo«  (/3j  -  a) 


A  V  sin  a  eat  $7 

vi  '-—rs - » -  mVa 


eat  (fit  -  a) 


(30) 


A 


*2  ~  Pi 
YP  i 


A 


2  • 


where  and  Pj  are,  respectively,  the  wave  angle  for  and  the  pressure 
behind  the  oblique  shock  OB . 

On  CD 


?!  •  A  -  0 


(31) 


On  DE 


Vy  m  -  AV  tin2a 

»>2  -  A  V  tin  a  cas  a  (32) 

A tin  a 
°1 


On  AE 


Vy  tin  a  -  Vj  eat  a  «  0 


£L  m  0 
dif 


(33) 


On  segnent  BC  a  difficulty  arises  in  that  the  details  of  the  shock  inter¬ 
action  are  not  known.  It  is  surmised  that  the  shock  segment  BC  is  a  smooth 
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curve,  matching  slopes  with  OB  at  B  and  CC“  atC  .  A  numerical  scheme  in 
which  the  shock  shape  is  determined  simultaneously  with  the  solution  of 
the  flow  equations,  in  the  manner  of  current  methods  of  solving  the  blunt 
body  problem,  would  seem  to  be  appropriate. 

If  the  wedge  angle  a  is  assuned  small,  shocks  OB  and  CC '  will,  in  the 
limit,  be  tangent  to  the  circle,  giving  the  geometry  of  Fig.  4.  In  this 
case  the  interaction  occurs  "outside"  the  unit  circle  and  the  conditions 
on  BC  are  a  simple  superpostion  of  the  conditions  on  AC  and  onBD.  That 
is,  A-Aj-Aj-l/y  and  v\mVi\  onCB. 


FIGURE  4 


7.  TOE  PRESSURE  DISTRIBUTION 


The  solution  for  a  may  now  be  written 


X  -  A,  — —  Rl 


|  log  (izillUi 
®,  *  V I  +  Ik  /  2 


- 0/7+  Aj)|l  teg  +  log  +(*c-«>jj. 


(34) 
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■where 


x  m  a*  tc  -  k*‘^C“ a* 


(35) 


This  solution  applies  inside  the  circle  ADE  .  Pig.  5  shows  a  plot  of  A 
versus  distance  along  wedge  surface  for  different  initial  Mach  nunbers. 


A 


8.  THE  VELOCITY  POTENTIAL 

For  the  small  angle  case,  the  entropy  gradient  is  of  the  order  «*  AV  and 
can  be  ignored.  The  coefficients  of  Eq.  (24)  can  then  be  formally  evalu¬ 
ated  .  We  define 
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tjf  m  fi  —  a 


and  write 


<f>  •  2  f  (n,r)  (C„  co»  ni/f  +  Dn  sin  nt/r) 


(36) 


Since  by  Eq.  (33)  -£JL_  -  0  at  ^  ®  and  tfr  -  w  ,  we  require  that  ^  be 
d\jt 

an  even  function  of  «|» ,  giving 


D„-0 


Eqs.  (30),  (3l)j  and  (32)  can  be  expressed  in  terms  of  i|r  as 


AV  tin  a  cot  \Ji 


dr 


■-  1 


L  v2\  cos  ifr  +  v32  sin  ^ 


0<*<f 


£<+<+ B 


(38) 


where  only  lowest  order  terms  have  been  retained. 


Formally 


i£ 

dr 


and  since  (Ref.  5) 


r'(f  T1'  1 ')  - ^  F (?r  "*'J ' ) 
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we  have 


f“  (l,f)  -  2  n2C„  eo«  n  ^  (39) 

Prcm  Eq.  (39),  using  standard  methods, 

tin  (n+1)  tin  (n-1) 
n  +  1  n - 1 

(*0> 

2(-l)n  cot  (n  +  1)  cot  ( n-1 )  1 

nJ-  1  n+l  n-1 


c„  — 


(-1)2  [1  +  (-1)"1 
trn2(n2  -  1) 


"21 

AVtina  - 

wn* 


9.  CONCLUSION 

Hie  pressure  distribution  and  the  potential  function  have  been  determined 
for  the  case  of  very  small  wedge  angle.  The  solution  for  larger  wedge 
angles  could  not  be  completely  obtained  owing  to  inability  to  specify  the 
boundary  values  in  the  shock  interaction  region.  However,  the  pressures 
to  either  side  of  the  unit  circle  are  constant  and  can  be  determined  from 
ordinary  shock  considerations.  It  is  felt  that  the  nature  of  the  pressure 
variation  on  the  wedge  surface  inside  of  the  unit  circle  is  not  markedly 
different  for  large  angles  from  that  given  by  the  small  angle  solution. 

Die  transient  foredrag  of  a  simple  wedge  wing  can  be  estimated  roughly  by 
noting  that  for  the  extreme  case  of  MaB=  1.25  the  pressure  coefficient  to 
the  rigit  of  E  is  approximately  six  percent  greater  than  that  to  the  left 
of  A .  The  mean  value  of  the  transient  foredrag  coefficient  would  then  be 
roughly  half  this  value,  that  is,  three  percent  of  the  final  steady  state 
foredrag  coefficient.  For  the  higher  Mach  nunbers  this  value  decreases, 
becoming  slightly  negative  for  M„.=  1.93  • 
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